Fluctuations of self-flattening surfaces 
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We study the scaling properties of self-flattening surfaces under global suppression on surface fluc- 
tuations. Evolution of self-flattening surfaces is described by restricted solid-on-solid type monomer 
deposition-evaporation model with reduced deposition (evaporation) at the globally highest (low- 
est) site. We flnd numerically that equilibrium surface fluctuations are anomalous with roughness 
exponent a ~ 1/3 and dynamic exponent zw — 3/2 in one dimension (ID) and q = (log) and 
zw — 5/2 in 2D. Stationary roughness can be understood analytically by relating our model to 
the static self-attracting random walk model and the dissociative dimer type deposition-evaporation 
model. In case of nonequilibrium growing/eroding surfaces, self-flattening dynamics turns out to be 
irrelevant and the normal Kardar-Parisi-Zhang universality is recovered in all dimensions. 
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Structural properties for fluctuating surfaces under 
thermal noise have been studied extensively j^. Equi- 
librium surfaces with proper surface tension are always 
rough in one dimension (ID) and display a roughening 
transition in two dimensions (2D) j^. Higher dimen- 
sional surfaces are always smooth. Surface roughness is 
well documented and classified as the Edwards- Wilkinson 
(EW) universality class ||^ . The E W class is generic and 
robust for equilibrium surfaces with local surface ten- 
sion. Only specific nonlinear contributions in nonequilib- 
rium growth processes may become relevant and drive the 
system into other universality classes, e.g., the Kardar- 
Parisi-Zhang (KPZ) universality class . 

In this paper, we introduce a new global mechanism 
to suppress surface fiuctuations, besides ordinary local 
surface tension. We call it self-flattening mechanism 
to reduce growing (eroding) probability at the globally 
highest (lowest) point on the surface. This global type 
suppression makes the surface less rough, which may 
bring forth new universality classes for equilibrium and 
nonequilibirum surfaces. Inclusion of suppression at all 
local extremal points leads to less interesting layer-by- 
layer growth processes and the steady state surfaces are 
always smooth with finite fluctuation width. 

We describe surface configurations in terms of inte- 
ger height variables {h{r)} at site r on a _D-dimensional 
hypercubic lattice. They are subject to the restricted 
solid-on-solid (RSOS) constraint, h{f-\-ei) — h{f) = 0, ±1 
with ii a primitive lattice vector in the i-th direction 
{i = 1, ...,D). The RSOS constraint effectively generates 
local surface tension which prevents indefinite growth of 
surface fiuctuations for finite systems. 

Evolution rule for the ordinary RSOS type monomer 
deposition-evaporation model is given as follows. First, 
select a site f randomly. Next, deposit a particle, 
h{r) — > h{r) + 1, with probability p or evaporate a parti- 
cle, h^r) — > h{r) — 1, with probability q = 1—p. Any de- 



position/evaporation attempt is rejected if it would result 
in violating the RSOS constraint. Equilibrium surfaces 
at p — q belongs to the EW class, while nonequilibrium 
growing/eroding surfaces at p q the KPZ class 

For self-flattening surfaces, we need a slight variation 
of the evolution rule to incorporate the global suppres- 
sion : only when deposition (evaporation) is attempted 
at the globally highest (lowest) site, the attempt is ac- 
cepted with probability u and rejected with probability 
1 — u. At M = 1, the ordinary RSOS model is recovered. 
The u = case is trivial, because the surface is conflned 
within initial surface width. 

We perform numerical simulations, starting from a flat 
surface of linear size L with periodic boundary condi- 
tions. We measure the surface fluctuation width W as 



(1) 



where (• • •) represents the ensemble average with equal 
weights. Therefore, our simulations at p = q correspond 
to the infinite temperature limit of equilibrium RSOS 
surfaces. The surface width satisfies the dynamic scaling 
relation 



W{L,t)^L''f{t/L' 



(2) 



where the scaling function f{x) — > const, for x ^ 1 and 
fix) r^xl^ {(3 = a/zw) for X < 1 

First, we report the numerical results for equilibrium 
surfaces (p — q). For ID, we run simulations for L = 
2^, ...,2^^ at u — 0.1, 0.3, 0.6, and 0.8, and average over 
at least 300 independent samples. In early time regime 
{t <C i^"'), the surface width grows with time, W ^ , 
and saturates to a finite value which increases with size. 
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In Fig. 1, we show the plot of In W versus \nt atu = 0.6 
for various system sizes. The growth exponent /3 is es- 
timated by a simple straight line fitting of early time 
data for the largest system size L = 2^^. Our estimate 
is /3 = 0.22(1) ~ 2/9. We also analyze the data at other 
values of u and find that (3 does not vary with u. 

In order to extract the stationary property, we average 
over data in the saturated regime {t ^ L^^) for given 
L to estimate Ws{L). For efficient estimation of a, we 
introduce effective exponents 



ae//(L) = ln [Ws{mL)/Ws{L)]/\nm, 



(3) 



where m is an arbitrary constant (here, we set m = 2). 

Effective exponents at various values of u are plotted 
in Fig. 2. Close to u = 1, our data show large correc- 
tions to scaling as expected, due to the presence of the 
EW fixed point [a — 1/2, P — 1/4) at w = 1. However, 
the asymptotic estimates seem to be independent of u. 
We estimate that a = 0.33(1) ~ 1/3 for all u. We check 
the dynamic scaling relation directly by plotting W/L" 
versus t/L^" in the inset of Fig. 1. Our data collapse 
very well with a = 1/3 and zw — 3/2 for all u, which 
are consistent with the above results. 

This set of scaling exponents form a new universality 
class, distinct from the EW and any previously known 
growth- type universality class. It implies that the self- 
flattening dynamics is a relevant perturbation to the EW 
fixed point in ID. Therefore, the continuum equation to 
describe self-flattening surfaces must contain a global- 
type nonlinear term. Further study in this direction is 
left for future research. 

In case of 2D EW surfaces, it is well known that the 
surface width grows logarithmically with time and its sat- 
urated value also increases logarithmically with size 
Especially, the saturated width Ws scales for large L as 
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where Kq is the effective coupling constant of the Gaus- 
sian model where equilibrium surface models flow into by 
renormalization group transformations The ordi- 

nary RSOS model at the infinite temperature (our model 
at u = 0) is known to take Kq ^ ^ 0.916 @J|. 
Assume the dynamic scaling relation similar to Eq. (|^) 

as 



W^{L,t) 
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where the scaling function g{x) — > const, for x ^ 1 and 
g{x) ~ x^/^"' for X <^ 1. Then, in early time regime 
{t -C i^"' ) , the surface width grows as 
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The amplitude ratio in Eqs.(|^) and (^) yields a value of 
the dynamic exponent zw. The EW surfaces take z^ = 2 
in all dimensions. 

We run simulations on LxL lattices with L = 2'^, 2'' 
at u = 0.1 and 0.5 and average over at least 300 inde- 
pendent samples. In Fig. 3(a), we plot against Int 
at u = 0.5. It shows a nice linear behavior in the early 
time regime. In Fig. 3(b), we plot against InL, which 
also shows a very nice linear behavior. We measure its 
slope and find that Kg — 0.92(1) for all u, which is very 
close to Kq. In contrast to the ID surfaces, the global 
suppression does not seem to change the asymptotic be- 
havior of the stationary surface roughness. As can be 
seen in Fig. 3, it seems to shift Wg only by a constant. 

We measure the amplitude ratio by comparing two 
slopes in Figs. 3(a) and (b). We estimate zw — 2.5(1) ~ 
5/2 for all u, which is clearly distinct from the EW value 
of 2. We also check the dynamic scaling relation of Eq.(||) 
by plotting VF^ — versus t/L^"^ in Fig. 3(a). Our 
data collapse reasonably well with zw = 5/2 for all u. 
Together with our ID results, we conclude that the self- 
flattening surfaces display a new type of scaling behavior 
and form a novel universality class. 

The partition function for equilibrium self-flattening 
surfaces can be written as 



Z = 



RSOS conf. 



(7) 



where the summation is over all height configurations sat- 
isfying the RSOS condition, (5 a temperaturelike param- 
eter, and ft-max (^min) thc globally maximum (minimum) 
height for a given configuration. 

Global suppression for self-flattening dynamics is sim- 
ply Metropolis type evolution algorithm with this parti- 
tion function to reach the equilibrium. Deposition (ero- 
sion) at the globally highest (lowest) site increases the en- 
ergylike term /imax — ^min by one unit and these attempts 
are accepted with Boltzmann type probability e^^ . Any 
other deposition (erosion) attempts are always accepted, 
because they do not increase the energylike term. Of 
course, all attempts resulting in violation of the RSOS 
constraint are rejected. By identifying u — e^^, our 
model for self-flattening surfaces is exactly the same as 
Metropolis evolution with the above partition function. 

Stationary property of this system can be understood 
analytically. In ID, this system is equivalent to the so- 
called static self-attracting (timid) random walks The 
surface can be mapped to the time trajectory of a random 
walker by identifying the height h(x) at site x with the 
walker position after x steps. The system size L becomes 
thc total number of steps and the RSOS constraint limits 
one-step hopping distance to or ±1. 

In ID, the energylike term is simply the number of dis- 
tinct sites visited by the random walker up to L steps. 
Random walk configurations with less visited sites are 
preferred. Such a random walker tends to visit previously 
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visited sites, so the walk is self-attractive. Its typical 
displacements are known rigorously to scale as 1,1/(^+2) 
P,p^ under the assumption that the visited sites form 
a compact cluster. In ID, the cluster is obviously com- 
pact, so the roughness exponent in our model should be 
a = 1/3 in ID. 

In 2D, the self-flattening surfaces are completely differ- 
ent from the self-attracting walks. The former deals with 
membrane fluctuations, while the latter polymer fluctu- 
ations. In order to understand the scaling behavior of 
the self-flattening surfaces, we investigate the intricate 
relation between our model and the dissociative dimer 
deposition-evaporation model in equilibrium 

In the dimer model, we deposit or evaporate particles 
only in a dimer form aligned along the surface. There 
is a global evenness conservation law that the number 
of particles at each height level must be conserved mod- 
ulo 2 This leads to a Boltzmann type factor in the 
partition function as 

RSOS conf. h 

where the product is over all possible height levels and 
Vh the number of particles at height level h. The dimer 
model corresponds to the z = — 1 case where only con- 
figurations obeying the evenness conservation law (all Vh 
are even) survive in the partition function. At z = 1, the 
model reduces to the ordinary monomer model. 

The self-flattening surfaces correspond to the z — 
limit. Each term inside the product picks up a factor of 
^ if Vh ^ 0, otherwise a factor of unity. The number 
of height levels with nonzero Vh (at least one particle) is 
^max ~ ^min- Therefore, the z — Q case is equivalent to 
the self-flattening surfaces at /3 = ln2. In fact, the Q-mer 
generalization corresponds to the /3 = In Q case |ll|,|l^ . 

/^From the Gaussian model type renormalization group 
argument, one can show that the 2D surface roughness 
is always logarithmic in the dimer model for —l<z<l 
(see Eq.(^) and its amplitude remains unchanged d). 
Our numerical results for all u are consistent with this. 
The dimer characteristics show up only in a form of cor- 
rections to scaling. Recently, it is suggested that the 
corrections to scaling should scale as In(lnL), which is 
confirmed for the dimer model at z = — 1 B. We find 
no evidence of this type of corrections to scaling in our 
model (z — 0) and the leading corrections are constants. 
The origin of this discrepancy between the z — Q and 
z = — 1 case is not fully understood as yet. 

Next, we consider noncquilibrium growing/eroding 
surfaces {p ^ q)- We run simulations for L ~ 2^, 2^^ 
for ID and L = 2^, ...,2^ for 2D at p = 1 with it = 0.5 and 
u — \ (ordinary RSOS). In Fig. 4, we plot InWs against 
InL and, in the inset, InM^ against Int for the largest 
system size in ID and 2D, respectively. We do not find 
any noticeable change of W ascribed to the global sup- 
pression. We estimate that a ~ 0.50(1) and /3 ~ 0.32(1) 



for ID and a ~ 0.40(1) and (3 ~ 0.24(1) for 2D, which are 
consistent with the results for the ordinary RSOS model 
Q . We conclude that the global suppression is irrelevant 
to nonequilibrium growing/eroding surfaces. 

In summary, we studied the scaling properties of the 
self-flattening surfaces in ID and 2D. Equilibrium sur- 
faces display dynamic scaling behavior distinct from the 
EW class and form a new universality class. Wc show 
that stationary roughness can be understood through 
mapping our model to self-attracting random walks in 
ID and dissociative dimer type deposition-evaporation 
model in 2D. In higher dimensions, the surfaces are al- 
ways smooth. In contrast, nonequilibrium self-flattening 
surfaces belong to the ordinary KPZ universality class. 
This implies that the self-flattening dynamics is strong 
enough to dominate over the EW type local surface ten- 
sion term, but weaker than the KPZ type nonlinear term. 
It would be very interesting to flnd a continuum-type 
equation to govern the self-flattening dynamics. 
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FIG. 3. (a) Plots of against Int at u = 0.5 for 2D self-fiattening equilibrium surfaces. The slope of the straight line yields 
the value of Kazw = 2.3(1). The inset shows the data collapse with Kg ~ 0.92 and zw = 5/2. (b) Plots of Ws^ against InL 
at « = 1 (ordinary RSOS), u = 0.5, and u = 0.1. The slopes of three straight lines yield the same value of Kg = 0.92(1) ~ Kq. 
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FIG. 4. Plots of In Ws against In L in ID and 2D nonequilibrium growii 
in W between the m = 1 (ordinary RSOS) and u = 0.5 (self-flattening) i 
for system sizes L = 2^^ (ID) and L = 2' (2D). Straight line fits yield 
/? = 0.24(1) for 2D. 
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